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Introduction 



The strip projection method is the most important way to generate quasiperiodic pat- 
terns with predefined local structure. We have obtained a very efficient algorithm for this 
method which allows one to use it in superspaces of very high dimension. A version of this 
algorithm for two-dimensional clusters and an application to decagonal two-shell clusters 
(strip projection in a 10-dimensional superspace) has been presented in |math-ph/ 0504036 
The program in FORTRAN 90 used in this case is very fast ( 700-800 points are obtained 
in 3 minutes ). 

We present an application of our algorithm to three-dimensional clusters. The phys- 
ical three-dimensional space is embedded into a 31-dimensional superspace and the strip 
projection method is used in order to generate a quasiperiodic packing of interpenetrating 
translated copies of a three-shell icosahedral cluster formed by the 12 vertices of a regular 
icosahedron (the first shell), the 20 vertices of a regular dodecahedron (the second shell) and 
the 30 vertices of an icosidodecahedron (the third shell). 

On a personal computer Pentium 4 with Fortran PowerStation version 4.0 (Microsoft 
Developer Studio) we obtain 400-500 points in 10 minutes. 

More details, bibliography and samples can be found on the website: 

http : / /fp cm5 . fizica.unibuc. ro~ncotf as / 



Computer program in FORTRAN 90 and MATHEMATICA 

! QUASIPERIODIC PACKINGS OF THREE-SHELL ICOSAHEDRAL CLUSTERS 
! ***** (ICOSAHEDRON + DODECAHEDRON + ICOSIDODECAHEDRON)***** 

! PLEASE INDICATE THE NUMBER OF POINTS YOU WANT TO ANALYSE 
INTEGER, PARAMETER : : N = 10000 

! PLEASE INDICATE THE DIMENSION M OF THE SUPERSPACE 
INTEGER, PARAMETER : : M = 31 
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INTEGER I, J, K, L, II, 12, 13, 14, JJ, JP, JPP 

REAL T, Rl, R2, R3, Dl, D2, D3, D4, AA 
REAL, DIMENSION (M) : : V, W, TR 
REAL, DIMENSION (3, 3) :: C5 
REAL, DIMENSION (3, M) :: B 

REAL, DIMENSI0N(1:M-3,2:M-2,3:M-1,4:M) :: S 
REAL, DIMENSION (N,M) :: P 
REAL, DIMENSION (N) :: X, Y, Z 

! PLEASE INDICATE THE RADIUS OF THE FIRST SHELL (ICOSAHEDRON) 
Rl = 1.0 

! PLEASE INDICATE THE RADIUS OF THE SECOND SHELL (DODECAHEDRON) 
R2 = 1.2 

! PLEASE INDICATE THE RADIUS OF THE THIRD SHELL (ICOSIDODECAHEDRON) 
R3 = 1.5 

! PLEASE INDICATE THE TRANSLATION OF THE STRIP YOU WANT TO USE 
TR = 0.1 

T = (l+SQRT(5.0))/2.0 
C5(l,l) = (T-D/2.0 
C5(l,2) = -T/2.0 
C5(l,3) = 1/2.0 
C5(2,l) = T/2.0 
C5(2,2) = 1/2.0 
C5(2,3) = (T-D/2.0 
C5(3,l) = -1/2.0 
C5(3,2) = (T-D/2.0 
C5(3,3) = T/2.0 
B = 0.0 

B(l,l) = Rl / SQRT(T+2.0) 

B(2,l) = T * Rl / SQRT(T+2.0) 

DO I = 2, 5 
DO J = 1, 3 
B(J,I) = SUM( C5(J,:) * B(:,I-1)) 
END DO 

END DO 

B(2,6) = Rl / SQRT(T+2.0) 
B(3,6) = T * Rl / SQRT(T+2.0) 

B(l,7) = R2 / SQRT(3.0) 
B(2,7) = R2 / SQRT(3.0) 
B(3,7) = R2 / SQRT(3.0) 
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DO I = 
DO J 

B(J,I) = 

END DO 
END DO 
B(l,12) 
B(2,12) 
BC3.12) 
DO I = 
DO J 

B(J,I) = 

END DO 
END DO 
B(l,17) 
DO I = 
DO J 

B(J,I) = 

END DO 
END DO 
B(2,22) 
DO I = 
DO J 

B(J,I) = 

END DO 
END DO 
B(3,27) 
DO I = 
DO J 

B(J,I) = 

END DO 
END DO 
PRINT* , 
PRINT* , 
PRINT* , 
PRINT* , 
PRINT* , 
PRINT* , 
PRINT* , 
DO J = 
PRINT* , 
END DO 
PRINT* , 



8, 11 
= 1, 3 

SUM( C5(J, :) * B(: ,1-1)) 



= R2 / SQRT(3.0) 

= -R2 / SQRT(3.0) 

= R2 / SORT (3.0) 
13, 16 
= 1, 3 

SUM( C5(J, :) * B(: ,1-1)) 



= R3 
18, 21 
= 1, 3 

SUM( C5(J, :) * B(: ,1-1)) 



= R3 
23, 26 
= 1, 3 

SUM( C5(J, :) * B(: ,1-1)) 



= R3 
28, 31 
= 1, 3 

SUM( C5(J, :) * B(: ,1-1)) 



'RADIUS OF THE FIRST SHELL (ICOSAHEDRON) IS ', Rl 

'RADIUS OF THE SECOND SHELL (DODECAHEDRON) IS ', R2 

'RADIUS OF THE THIRD SHELL (ICOSIDODECAHEDRON) IS ', R3 

'STRIP TRANSLATED IN SUPERSPACE WITH THE VECTOR OF COORDINATES:' 

TR 

'COORDINATES OF THE POINTS OF THREE-SHELL ICOSAHEDRAL CLUSTER' 
' (UP TO A SYMMETRY WITH RESPECT TO THE ORIGIN) : ' 
1, M 

J, B(1,J), B(2,J), B(3,J) 

'PLEASE WAIT A FEW MINUTES OR MORE,& 
DEPENDING ON THE NUMBER OF ANALYSED POINTS' 



S = 
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DO II = 1, M-3 
DO 12 = 11+1, M-2 
DO 13 = 12+1, M-l 
DO 14 = 13+1, M 
DO Dl = -0.5, 0.5 
DO D2 = -0.5, 0.5 
DO D3 = -0.5, 0.5 
DO D4 = -0.5, 0.5 

AA = Dl * ( B(1,I2) * B(2,I3) * B(3,I4) + & 
B(2,I2) * B(3,I3) * B(1,I4) + & 
B(3,I2) * B(1,I3) * B(2,I4) - & 
B(3,I2) * B(2,I3) * B(1,I4) - & 

B(1,I2) * B(3,I3) * B(2,I4) - & 
B(2,I2) * B(1,I3) * B(3,I4) ) - & 
D2 * ( B(1,I1) * B(2,I3) * B(3,I4) + & 

B(2,I1) * B(3,I3) * B(1,I4) + & 
B(3,I1) * B(1,I3) * B(2,I4) - & 
B(3,I1) * B(2,I3) * B(1,I4) - & 

B(l.Il) * B(3,I3) * B(2,I4) - & 
B(2,I1) * B(1,I3) * B(3,I4) ) + & 
D3 * ( B(1,I1) * B(2,I2) * B(3,I4) + & 

B(2,I1) * B(3,I2) * B(1,I4) + & 
B(3,I1) * B(1,I2) * B(2,I4) - & 
B(3,I1) * B(2,I2) * B(1,I4) - & 

B(1,I1) * B(3,I2) * B(2,I4) - & 
B(2,I1) * B(1,I2) * B(3,I4) ) - & 
D4 * ( B(1,I1) * B(2,I2) * B(3,I3) + & 

B(2,I1) * B(3,I2) * B(1,I3) + & 
B(3,I1) * B(1,I2) * B(2,I3) - & 
B(3,I1) * B(2,I2) * B(1,I3) - & 

B(l.Il) * B(3,I2) * B(2,I3) - & 
B(2,I1) * B(1,I2) * B(3,I3) ) 

IF ( AA > S(I1, 12, 13, 14) ) S(I1 , 12, 13, 14) = AA 
END DO 

END DO 
END DO 
END DO 

IF( S(I1, 12, 13, 14) == 0) S(I1,I2,I3,I4) = N * SUM( B(l,:) ** 2 ) 
END DO 

END DO 
END DO 
END DO 
P = 

P( 1, :) = ANINT( TR ) 
K = 1 
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L = 

JP = 

DO I = 1, N 

IF( I <= K ) 

V = P(I, : ) 

JJ = 1 

JPP = 

DO II = 1, 



THEN 
- TR 



M-3 



DO 12 = 11+1, M-2 
DO 13 = 12+1, M-l 
DO 14 = 13+1, M 
AA = V(I1) * ( B(1,I2) * B(2,I3) * B(3,I4) + & 
B(2,I2) * B(3,I3) * B(1,I4) + & 
B(3,I2) * B(1,I3) * B(2,I4) - & 
B(3,I2) * B(2,I3) * B(1,I4) - & 

B(1,I2) * B(3,I3) * B(2,I4) - & 

* B(1,I3) * B(3,I4) ) - & 
B(1,I1) * B(2,I3) * B(3,I4) + & 

B(2,I1) * B(3,I3) * B(1,I4) + & 

* B(1,I3) * B(2,I4) - & 

* B(2,I3) * B(1,I4) - & 
B(l.Il) * B(3,I3) * B(2,I4) - & 

* B(1,I3) * B(3,I4) ) + & 
B(1,I1) * B(2,I2) * B(3,I4) + & 

B(2,I1) * B(3,I2) * B(1,I4) + & 

* B(1,I2) * B(2,I4) - & 

* B(2,I2) * B(1,I4) - & 
B(l.Il) * B(3,I2) * B(2,I4) - & 

* B(1,I2) * B(3,I4) ) - & 
B(1,I1) * B(2,I2) * B(3,I3) + & 

B(2,I1) * B(3,I2) * B(1,I3) + & 
B(2,I3) - & 
B(1,I3) - & 

* B(3,I2) * B(2,I3) - & 



B(2,I2 
V(I2) * 

B(3,I1 
B(3,I1 

B(2,I1 
V(I3) * 

B(3,I1 
B(3,I1 

B(2,I1 
V(I4) * 

B(3,I1 
B(3,I1 



B(2,I1) * 



B(1,I2) * 
B(2,I2) * 
B(1,I1) 
B(1,I2) * 



( AA < -S(I1,I2,I3,I4) 



IF 
IF 
END DO 



( AA == -S(I1, 12, 13, 14) 



END DO 
END DO 
END DO 



B(3,I3) ) 

.OR. AA > S(I1, 12, 13, 14) ) JJ = 
.OR. AA == S(I1, 12, 13, 14) ) JPP = 1 



IF ( J J .EQ. 1 ) THEN 
XP = SUM( V * B(l, :) ) 

YP = SUM( V * B(2, :) ) 
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ZP = SUM( V * B(3, :) ) 

13 = 1 
IF ( L > 0) THEN 

DO J = 1, L 

IF( XP == X(J) .AND. YP == Y(J) .AND. ZP == Z(J) ) 13 = 
END DO 

ELSE 
END IF 

IF( 13 == 1) THEN 
IF( JPP .EQ. 1 ) JP = JP + 1 
L = L + 1 
X(L) = XP 

Y(L) = YP 
Z(L) = ZP 
ELSE 

END IF 
DO II = 1, M 
DO 12 = -1, 1 
W = P(I, :) 
W(I1) = W(I1) + 12 
13 = 
DO J = 1, K 

IF( ALL(W .EQ. P(J, :)) ) 13 = 1 
END DO 

IF( 13 == .AND. K < N ) THEN 
K = K + 1 
P(K, : ) = W 
ELSE 
END IF 

END DO 
END DO 
ELSE 

END IF 
ELSE 
END IF 
END DO 

PRINT*, ' NUMBER OF ANALYSED POINTS : K 
PRINT*, ' NUMBER OF OBTAINED POINTS : ', L 

PRINT*, ' NUMBER OF POINTS LYING ON THE FRONTIERE OF THE & 

STRIP:', JP 

PRINT*, 'PLEASE INDICATE THE NAME OF A FILE FOR RESULTS' 
WRITE (4, 98) 

98 FORMAT ( ' Show [Graphics3D [{ PointSize [0 . 01] , { ' ) 
DO J = 1, L-l 

WRITE(4,99) X(J), Y(J) , Z(J) 
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99 FORMAT ( ' Point [{ 'FIO . 5 ' , ' F10 . 5 ' , ' FIO . 5 , ' }] , ') 
END DO 

WRITE(4,100) X(L), Y(L) , Z(L) 

100 FORMAT ( 'Point [{'FIO. 5' , 'FIO. 5' , 'FIO. 5, '}] ') 
WRITE(4,101) 

101 FORMAT ( ' }} ]] ') 

PRINT*, '* OPEN THE FILE CONTAINING THE RESULTS WITH & 

"NotePad" ' 

PRINT*, '* SELECT THE CONTENT OF THE FILE ("Select All") & 

AND COPY IT ("Copy") ' 
PRINT*, '* OPEN "MATHEMATICA", PASTE THE COPIED FILE, & 

AND EXECUTE IT ("Shif t+Enter") ' 

END 
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